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Abstract 

Background: Few studies have been conducted in Pakistan to determine the efficacy of chloroquine and 
sulphadoxine-pyrimethamine (SP), which remain in use as treatment for Plasmodium vivax and in combination with 
artesunate to treat Plasmodium falciparum, respectively. In this study, samples from several sites across Pakistan 
were characterized to determine prevalence of molecular resistance markers in the P. falciparum chloroquine 
resistance transporter ipfcrt), multidrug resistance ipfmdrl), dihydrofolate reductase ipfdhfr) and dihydropteroate 
synthase (pfdhps) genes and the origin of chloroquine-resistant P. falciparum parasites. 

Methods: Microscopy-confirmed malaria parasite-positive blood samples from 801 patients across the country were 
collected in 201 1. Of these, 171 infections were identified by polymerase chain reaction (PGR) as P. falciparum and 
analysed by pyrosequencing for mutations conferring chloroquine resistance {pfcrt codons 72-76), multidrug 
resistance ipfmdrl N86Y, Y184F, S1034C, N1042D and D1246Y), pyrimethamine resistance [pfdhfr, C50R, N51I, C59R, 
S108N and I164L) and sulphadoxine resistance [pfdhps, S436A, A437G, K540E, A581G and A613T/S). pfmdrl gene 
copy number variation was determined by real-time PGR, and microsatellites flanking the pfcrt locus were typed to 
determine the origin of the chloroquine-resistant haplotype. 

Results: The pfcrt K76T mutation was found in all samples as part of the S72/V73/M74/N75/T76 (SVMN"0 haplotype. 
Microsatellites flanking pfcrt showed high similarity to the signature found in India and Papua New Guinea, pfmdrl 
N86Y was found in 20% of samples and all samples harboured a single copy of the pfmdr] gene. The pfdhfr double 
mutation G59R + S108N was present in 87% of samples while the pfdhfr tuple mutant (N51 1 C59R + S108N) was 
not detected. Pfdhps A437G was found in 60% of samples. Pure pfdhps K540E was rare, at 4%, but mixed genotype 
540 K/E was found in 77% of samples. Similarly, pure pfdhps A581G was found in 4% of the isolates while mixed 
581 A/G was found in 39% of samples. 

Conclusions: These results suggest an emerging problem with multidrug resistant P. falciparum in Pakistan. The 
chloroquine resistance genotype has reached complete fixation in the population, with a microsatellite pattern 
indicative of a selective sweep. Moreover, the prevalence of mutations in both pfdhfr and pfdhps, albeit without the 
presence of the pfdhfr tuple mutant, indicates that continued monitoring is warranted to assess whether SP 
remains efficacious as a partner drug for artesunate for the treatment of P. falciparum. 
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Background 

In Pakistan, an estimated 500,000 episodes of malaria infec- 
tion occur annually [1]. Although the majority of these 
cases is attributed to Plasmodium vivax, studies indicate 
that Plasmodium falciparum in Pal<istan has been on the 
rise over the past few decades [2]. The proportion of mal- 
aria infections attributed to P, falciparum rose from 34 to 
54% between 1987 and 1990 in north-west Pakistan [3], 
and the frequency of P, falciparum among microscopy- 
positive cases rose from 45% in 1995 to 68% in 2006 in the 
provinces of Balochistan and Sindh [4,5]. In 2010, out of a 
total of 240,591 reported malaria cases in Pakistan, 73,857 
(31%) were P. falciparum [6]. 

The first-line treatment of uncomplicated malaria of 
undetermined species in Paldstan is chloroquine. Chloroquine- 
resistant P. falciparum was first reported in Pakistan in 
1984 [3] and subsequent in vivo studies have confirmed 
that chloroquine is not an efficacious treatment for 
P. falciparum in Central Asia. A survey conducted in 
western Pakistan in Afghan settlements during 1994-1995 
reported in vivo chloroquine resistance ranging from 18 
to 62% [7]. Both standard and extended-dose courses of 
chloroquine resulted in high levels of treatment failure in 
clinical studies conducted in Afghani refugees in Pakistan 
[8]. Other clinical studies conducted in refugee communi- 
ties throughout the country indicated that chloroquine re- 
sistance increased as much as 15% in exposed populations 
in a single year, rising five- to six-fold between 1982 and 
1992 [2]. 

In response to the declining efficacy of chloroquine, 
the World Health Organization (WHO) recommended 
artesunate plus sulphadoxine-pyrimethamine (AS + SP) 
as the first-line choice of treatment for uncomplicated P. 
falciparum malaria in Pakistan [9]. However, a 56% 
treatment failure rate with SP monotherapy was 
reported in Balochistan province in a study conducted 
between 2001 and 2005 [10], indicating that SP may be 
compromised as an effective partner drug. 

Molecular surveillance of drug resistance has been 
used extensively to monitor sensitivity to chloroquine 
and SP. A point mutation causing an amino acid change 
from K to T at codon 76 of the P, falciparum chloro- 
quine resistance transporter gene {pfcrt) confers resist- 
ance to chloroquine in vitro [11,12] and is strongly 
associated with treatment failure in vivo [13]. The role 
of these mutations in chloroquine resistance has also 
been confirmed by transfection studies [11]. A small 
number of independent pfcrt resistant lineages has been 
found: two in South America, one in Papua New Guinea 
[14] which has also been found in Asia [15,16], one ob- 
served in the Philippines, and one in Southeast Asia 
which spread to Africa [14,17]. 

Mutations in the dihydropteroate synthase (pfdhps) and 
dihydrofolate reductase (pfdhfr) genes of P, falciparum 



confer resistance to SP. Point mutations at positions N51I, 
C59R, S108N and I164L in pfdhfr are associated with pyri- 
methamine resistance in vitro [18-20] and treatment fail- 
ure in vivo [21-23]. Resistance to sulfas is conferred by 
mutations at codons A437G, K540E, A581G and A613T/S 
of the pfdhps gene, confirmed by in vitro [24-28] and 
in vivo [22-24] studies. Resistance to SP occurs in a step- 
wise fashion, with an increasing number of mutations in 
pfdhfr and pfdhps contributing to an increased risk of 
treatment failure [29]. Infections carrying the pfdhfr triple 
mutant (S108N + N51I + C59R) are significantly more 
likely to fail SP treatment than infections with fewer pfdhfr 
mutations [30], and the pfdhfr /pfdhps quintuple mutant 
(pfdhfr triple mutant S108N + N51I + C59R plus the pfdhps 
double mutant A437G + K540E) is highly predictive of 
clinical failure in Africa [23,29,31]. 

Elevated copy number of the pfmdrl gene is associated 
with susceptibility of P, falciparum to artemisinin, mef- 
loquine, and lumefantrine [32-34], and has been used as 
a surveillance tool for artesunate-mefloquine resistance 
in Southeast Asia [33,35]. Mutations in the pfmdrl gene 
leading to the substitution of amino acids including 
N86Y, Y184F, S1034C, N1042D, and D1246Y also modulate 
in vitro susceptibility to a number of drugs [36], although 
their relevance in conferring clinical resistance remains 
uncertain. 

In Pakistan, previous molecular surveys of drug resist- 
ance have focused on one or two sites in a small number 
of geographic regions. In these surveys, the most recent 
of which was conducted in 2007, high levels of pfcrt 76T 
have been reported [1,15,37], and multiple mutations in 
the pfdhfr and pfdhps genes have been found [1,37,38]. 
In this study, the prevalence of molecular markers of 
anti-malarial resistance was investigated in a compre- 
hensive survey over 14 sites in four provinces of Pakistan 
and in the capital city of Islamabad. Collected in 2011, 
these samples provide information on the current distri- 
bution of drug-resistant genotypes in the pfcrt, pfmdrl, 
pfdhfr and pfdhps genes five years after the implementa- 
tion of artesunate + SP as the first-line treatment for un- 
complicated P, falciparum. 

Methods 

Study sites and ethics 

Government and private hospitals in 25 cities from four 
provinces with the highest burden of malaria representing 
all four provinces (Khyber Pakhtunkhwa province, Sindh 
province, Balochistan province and Punjab province) and 
a hospital in the capital city, Islamabad, were invited to 
participate in the molecular survey in 2011. Although the 
Federally Administered Tribal Areas have among the 
highest burdens of malaria in the country, they were 
excluded from this study because political instability 
and violence [39] make it difficult to estabUsh sample 
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of Medicine in Baltimore, MD, USA. DNA was extracted 
from the blood-impregnated filter papers using the 
QIAmp 96 DNA kit (Qiagen, Valencia, CA, USA). Plas- 
modium species {P. vivax and P, falciparum) were con- 
firmed by nested PGR amplification of the small subunit 
ribosomal ribonucleic acid (ssrRNA) genes by using the 
primers and thermal cycler conditions as previously de- 
scribed [41]. PGR products were visualized by 2-2.5% 
agarose gel electrophoresis using the Bio-Rad gel doc 
system. Extracted DNA was stored at -20°G for further 
analysis. 

Pyrosequencing of pfdhfr, pfdhps, pfcrt, and pfmdrl 
codons 

PGR for pyrosequencing of pfdhfr codons G50R, N51I, 
G59R, S108N and I164L and pfdhps codons S436A, 
A437G, K540E, A581G and A613T/S was conducted 
using previously described cycling conditions [42] with a 
few modifications. Nested PGR reactions were prepared 
in 25 (il PGR reaction volume, contained 1 \A DNA tem- 
plate, Ix PGR buffer (Qiagen, Valencia, GA, USA), 
0.2 mM dNTPs (Invitrogen) and 0.05 units HotStar Taq 
DNA polymerase (Qiagen, Valencia, GA, USA). The sub- 
stitution of HotStar Taq required an initial incubation 
step of at 95°G for 15 min for all PGR reactions. PGR 




Figure 1 2011 collection sites with malaria-positive isolates: 14 cities in four provinces and the capital of Pakistan. 

V ) 



collection. Of the facilities that were contacted, hospi- 
tals in 14 cities shared samples and were included in this 
study (Figure 1). Symptomatic patients of all age groups 
were enrolled and 3 mL of venous blood was drawn from 
each patient. Gonsent was obtained from patients or from 
children's parents or guardians. Information on patient 
age and gender was also collected. The study was ap- 
proved by the institutional review board of Quaid-i-Azam 
University, Islamabad, Pakistan. 

Microscopy and sample collection 

Thick and thin blood films of patients with suspected 
malaria were stained with 10% Giemsa solution and ex- 
amined at 1,000 X magnification under oil immersion by 
a laboratory technician or technologist trained in malaria 
diagnosis according to WHO guidelines [40]. Approxi- 
mately 50 (il of each blood sample from 801 malaria 
microscopy-positive and 30 microscopy-negative sam- 
ples was applied to Whatman 3MM filter paper. Filter 
papers were air-dried overnight and stored at room 
temperature in individual plastic bags with desiccant. 

DNA extraction and speciation 

Molecular analysis was conducted at the Howard 
Hughes Medical Institute/University of Maryland School 
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was performed busing a BioRad tetrad Thermal cycler 
(Bio-Rad, Hercules, CA). PGR primers and cycling con- 
ditions for pyrosequencing of pfcrt codons 72-76 and 
pfmdrl codons N86Y, Y184F, S1034C, N1042D and 
D1246Y are described on the investigators website [43]. 

Pyrosequencing for all genes was performed using 
PyroMark® Q96 MD pyrosequencer, using the protocol 
previously described [42]. Between 2 and 7 \A of second- 
ary PGR product or sequence-specific amplified positive 
control DNA provided by Malaria Research Reagent Re- 
source (MR4) (Manassas, Virginia, USAJ was added to 
each pyrosequencing reaction, depending on the on the 
PGR product yield. Single nucleotide polymorphisms 
(SNPs) were called using PyroMark® Q96 MD 
pyrosequencing software version 1.2 (Qiagen) in allele 
quantification mode (AQ) for all SNPs except for pfcrt 
72-76 and pfmdrl N86Y, which were analysed in se- 
quence analysis mode (SQA). Single peak signals of at 
least 30 RLU (relative luminescence units) were consid- 
ered suitable for allele quantification. Pyrosequencing 
was repeated with an adjusted amount of PGR product 
for samples failing to produce a SNP call because of too 
much DNA (saturating the pyrogram) or insufficient 
DNA (producing a weak signal). 

SNP calls obtained by pyrosequencing were adjusted using 
standard curve equations derived from pyrosequencing of 
mixtures of control DNA strains with known proportions 
for each allele. To account for machine error, the lower 
cutoff for inclusion of minor alleles in a mixed infec- 
tion was set at 10%. pfcrt K76T and pfmdrl N86Y could 
not be absolutely quantified using pyrosequencing so a 
standard curve could not be applied. The minor allele fre- 
quency cutoff for these two codons was set at 20% to ac- 
count for both machine error and lack of adjustment by 
standard curves. 

Isolates were considered to be pure' mutants if only 
the resistance-conferring codon was observed, and 
mixed' if both the wild-type and mutant allele were 
detected. Similarly, haplotypes (combinations of poly- 
morphisms) were considered pure' mutants if wild-type 
alleles were not detected in any of the SNP positions, 
and mixed if both wild-type and mutant alleles were 
detected in at least one SNP position. 

pfmdrl copy number variation 

pfmdrl copy number was assessed by TaqMan real-time 
PGR (ABI sequence detector 7700; Applied Biosystems, 
Warrington, UK) using previously described primers, 
probes, and reaction conditions [33]. Every qPGR run 
contained genomic DNA extracted from P. falciparum 
strains of 3D7, Dd2 and Kl with pfmdrl copy numbers, 
of 1, 2, and 3 respectively. Each sample in this assay was 
run in triplicate. Assays were repeated if one of the fol- 
lowing three results was obtained: AAGt spread > 1.5; Gt 



values > 35; or copy number estimates between 1.3 and 
1.6 [33]. 

Microsatellite genotyping 

Six microsatellites, located 2.8, 4.3, 10.8, and 29.3 kb up- 
stream of pfcrt and 0.6 and 10.4 downstream of pfcrt on 
chromosome 7, were used for genotyping. These 
microsatellites were amplified using heminested PGR 
and fluorescently labeled primers [44]. Specifics of PGR 
reactions, thermal cycling conditions, and fragment size 
calling were carried out as described by Laufer et al 
[45]. In addition to samples collected in this study, four 
P, falciparum DNA samples collected between 2005 and 
2007 in a previous study in Pakistan [15] and four col- 
lected in Dagua Province in Papua New Guinea in 2010 
were genotyped at all six microsatellite loci to assess 
shared ancestry of the pfcrt locus. 

Results 

A total of 801 microscopy-positive samples was collected 
from 14 sites. Of these isolates, 128 samples were identi- 
fied by PGR as P, falciparum and an additional 43 sam- 
ples contained both P, falciparum and P, vivax 
infections (Table 1) for a total of 171 samples containing 
P, falciparum. Sixty three percent (108 samples) were 
collected from male subjects and 37% (63 samples) were 
collected from females. Subject age ranged from three 
months to 75 years, with a median of 24 years. 

Pyrosequencing of pfdhfr and pfdhps genes harbouring 
mutations conferring drug resistance was carried out in 
171 PGR-positive P, falciparum isolates. The pfdhfr 
triple mutant N51I + G59R + S108N was detected in 7 
samples, all carrying a mixture of wild-type and mutant 
alleles (Table 2). The pfdhfr double mutant G59R + S108N 
was found in 148 samples (87%) (Table 2). The pfdhps 
mutant codon A437G was observed in 61 samples 
(51.5%), with few mixed infections, pfdhps K540E was 
found in six samples, and mixed K/E infections were 
found in 131 samples (77%). pfdhps A581G was observed 
in five samples while 66 samples (39%) had the mixed 
A/G genotype, pfdhfr I pfdhps G59R + S108N + A437G 
was found in 88 isolates (51%). 

The distribution of common pfdhfr and pfdhps haplo- 
types (combinations of polymorphisms) by region is 
shown in Table 3. Prevalence of pfdhfr double mutant 
G59R + S108N ranged from 72% in Punjab (in 31 of 43 
samples) to 95% in Khyber Pakhtunkhwa (in 87 of 91 
samples), pfdhfr C59R SlOm pfdhps A437G was 
present in 44-58% of samples collected in Punjab, Khyber 
Pakhtunkhwa, and Balochistan (sample sizes of 43, 91, 
and 33, respectively). Only two P. falciparum samples 
each was collected from Islamabad and Sindh, and all four 
carried pfdhfr G59R + S108N + pfdhps A437G. 
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Table 1 Plasmodium falciparum samples collected from 14 
sites In Pakistan 

Province City PCR-confirmed 

P. falciparum (N) 



Khyber Pakhtunkhwa 


Bannu 


39 




Hangu 


19 




Peshawar 


4 




Thai 


29 


Punjab 


Bhakkar 
Lahore 


1 
1 




Mianwali 


28 




Muzaffarghar 


13 




Rawalpindi 


0 


Balochistan 


Quetta 


26 




Zhob 


7 


Sindh 


Hyderabad 


0 




Karachi 


2 


Capital 


Islamabad 


2 


Total 




171 



carried the mixed Y/F codon at position Y184F of the 
pfmdrl gene and the remaining samples carried the 
wild-type Y codon. One hundred and sixty-five (96%) 
isolates carried the mixed N/D polymorphism at pos- 
ition N1042D. All samples were wild-type for codons 
S1034C and D1246Y. Real time PGR to determine copy 
number variation in pfmdrl gene revealed that all iso- 
lates harboured a single copy of the pfmdrl gene. 

All six microsatellite loci amplified successfully in 160 
of 171 P falciparum samples. Three isolates were mixed- 
strain infections with multiple alleles at each locus and 
were excluded from further analysis. Three haplotypes 
were observed in the remaining sample set of 157: the 
predominant haplotype, found in 137 samples, and two 
haplotypes in 12 and eight samples, respectively, each 
differing from the predominant haplotype at one micro- 
satellite locus (Table 4). These pfcrt flanking microsatel- 
lite haplotypes were identical or nearly identical 
(differing in 1 locus) to those found in isolates with the 
SVMNT haplotype from Pakistan [15] and Papua New 
Guinea (Table 4). 



All 171 samples analysed were amplified successfully 
for pfcrt. DNA pyrosequencing of the secondary PGR 
products from these samples confirmed that all of the 
amplified samples carried codon K76T (nucleotide se- 
quence AG A). All samples were of the SVMNT haplo- 
type for pfcrt 72-76 (Table 2). 

Pfmdrl N86Y pure mutants were identified in six (4%) 
of the tested isolates, 28 (16%) were identified as mixed 
N/Y infections (Table 2). Forty-three (25%) isolates 



Discussion 

In this study, polymorphisms in the pfdhfr, pfdhps, pfcrt 
and pfmdrl genes were investigated in P. falciparum 
samples collected throughout Pakistan to determine the 
current extent of resistance to chloroquine and 
antifolates, both of which are still used as part of first- 
line treatment of malaria in Pakistan. 

The complete saturation of pfcrt K76T in all parts of 
the country indicates that chloroquine-resistant P, fal- 
ciparum is fixed in Pakistan. These findings corroborate 



Table 2 Number (N) and prevalence (%) of pure/mixed pfcrt, pfmdrl , pfdhfr and pfdhps mutant alleles conferring 
resistance to chloroquine and sulphadoxine-pyrimethamine in Plasmodium falciparum isolates from Pakistan 







Pure 


Mixed 


Gene 


Allele or haplotype 


N (%) 


N (%) 


pfcrt 


SVMNT/ K76T 


171 (100) 


0(0) 


pfmdrl 


N86Y 


6(4) 


28 (16) 




Y184F 


0(0) 


43 (25) 




S1034C 


0(0) 


0(0) 




N1042D 


0(0) 


165 (96) 




D1246Y 


0(0) 


0(0) 


pfdhfr 


pfdhfr triple (N51 1 + C59R + SI 08N) 


0(0) 


7(4) 




pfdhfr double (N51I + S108N) 


11 (6) 


3(2) 




pfdhfr double (C59R + S108N) 


148 (87) 


9(5) 


pfdhps 


pfdhps double (A437G + K540E) 


4(2) 


77 (45) 




pfdhps A437G 


103 (60) 


3(2) 




pfdhps K540E 


6(4) 


131 (77) 




pfdhps A581G 


5(3) 


66 (39) 


pfdhfr + pfdhps 


pfdhfr + pfdhps (C59R + S108N + A437G) 


88 (51) 


6(4) 
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Table 3 Number and prevalence of SP resistance-associated haplotypes in pfdhfr and pfdhps in Pakistan, by province 

Number and prevalence (%) of pfdhfr, pfdhps, and pfdhfr + pfdhps alleles 
pfc//7fr N51I + C59R+ pfcf/ifr N51I + SI 08N pfdhfr C59R+ pfdhps M37G+ pfdhfr C59R + SWSN+ pfdhps 







S108N 




S108N 




K540E 


A437G 




n 


Mixed* 


Pure** 


Pure 


Mixed 


Pure 


Mixed 


Pure 


Mixed 


Balochistan 


33 


4(12) 


1 (3) 


26 (79) 


1 (3) 


2(6) 


16 (49) 


19 (58) 


0 




91 


0 


2(2) 


87 (95) 


1 (1) 


1 (1) 


32 (35) 


46 (51) 


1 (1) 


Islamabad 


2 


0 


0 


2 (100) 


0 


0 


2 (100) 


2(1) 


0 


Punjab 


43 


3(7) 


8 (19) 


31 (72) 


1 (2) 


1 (2) 


25 (58) 


19 (44) 


5(12) 


Sindh 


2 


0 


0 


2 (100) 


0 


0 


2 (100) 


2 (100) 


0 


Total 


171 


7(4) 


11 (6) 


148 (87) 


3(2) 


4(2) 


77 (45) 


88 (52) 


6(4) 



*Mixed pfdhfr N51I + S108N not detected. 
*^Pure pfdhfr triple mutant not detected. 
*^*KPK: Khyber Pakhtunkhwa. 



previous molecular studies reporting near fixation of 
chloroquine resistance [15,37] and clinical studies 
reporting high treatment failure of chloroquine [2,7,8]. 
Chloroquine resistance has also been documented in 
several neighbouring countries, including Iran, India, 
China and Afghanistan [46-50]. Although chloroquine is 
recommended for P. vivax infections only, P, falciparum 
infections are often exposed to this treatment as well. 
Presumptive diagnosis based on clinical grounds, im- 
proper diagnosis, lack of diagnostic facilities, and empir- 
ical treatment are common practices in resource-limited 
countries such as Pakistan [51,52], and these factors 



likely result in the continued mistreatment of P. falcip- 
arum infections with chloroquine. Such ineffective treat- 
ment of chloroquine-resistant P, falciparum with 
chloroquine can lead to increased morbidity and mortal- 
ity if infections do not respond to the drug. Use of 
chloroquine to treat P, falciparum has also been postu- 
lated as contributing to the rapid growth of falciparum 
malaria in Pakistan [2], as infections may not clear and 
continue to circulate in the population unchecked. 

Typing of microsatellites surrounding the pfcrt locus 
indicates that a selective sweep of chloroquine-resistant 
P, falciparum occurred in Pakistan, and that the 



Table 4 Microsatellite haplotypes of Plasmodium falciparum parasites In the region flanking the pfcrt locus on 
chromosome 7 



Microsatellites flanking pfcrt locus^ 



Country of origin 


Number or strain 


pfcrt haplotype 


-29 


-11 


-4 


-3 


+1 


+10 


Pakistan (current study) 


137 


SVMNT 


147 


168 


230 


178 


151 


203 




12 


SVMNT 








182 








8 


SVMNT 












197 


Pakistan** 


4 


SVMNT 














Papua New Guinea 


2 


SVMNT 
















2 


SVMNT 


149 












Papua New Guinea 


DIO 


CVMNK 




170 




166 




193 


Tliailand 


Kl 


CVIET 






228 




149 


200 


Papua New Guinea 


FC27 


CVIET 




170 




166 


149 


193 


Indocliina 


Dd2 


CVIET 




176 


228 




149 


200 


Vietnam 


VIS 


CVIET 




176 


228 




149 


200 


Indocliina 


W2 


CVIET 




176 


228 




149 


200 


Africa 


3D7 


CVMNK 


149 


180 




180 


149 




Brazil 


7G8 


CVMNT 


149 


174 






157 


200 


Honduras 


HB3 


CVMNK 


149 


182 




184 


147 


193 


Sierra Leone 


D6 


CVMNK 


145 


186 


232 


186 


149 


200 



^Numbers indicate distance in kilobases upstream (-) and downstream (+) of the pfcrt locus. 
**Reference [14]. 
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SVMNT allele in Pakistan shares its ancestry with para- 
sites in Papua New Guinea. Previous studies have shown 
that the SVMNT allele of the pfcrt locus appears to have 
swept through Papua New Guinea [14], multiple parts of 
India [16,46] and a small number of sites Pakistan [15]. 
Taken together with the findings in this study, these 
results suggest that the spread of SVMNT may have 
been a single sweep event through a large region 
encompassing parts of the South Pacific and South Asia. 
The possibility that such a large sweep could occur again 
in the context of artemisinin resistance in this region 
should be carefully considered. 

The current first-line treatment for confirmed P. fal- 
ciparum malaria in Pakistan is artesunate + SP. The 
pfdhfr C59R + S108N double mutant was the most fre- 
quent pyrimethamine-resistant haplotype found, as has 
been shown in previous studies in Pakistan [1,37], India 
[47,48], Iran [38,49], and Sri Lanka [50]. The majority of 
samples also carried a mutation at position A437G and a 
number of samples harboured mixed infections at highly 
resistant codons K540E and A581G in the pfdhps gene. 
Although the clinically significant pfdhfr triple mutation 
was rare (and present only in mixed wild-type/mutant 
infections) and pfdhfr/pfdhps quintuple mutants were 
not detected, the presence of multilocus resistance hap- 
lotypes and the presence of highly resistant sulphadoxine 
alleles indicate that clinically relevant SP resistance 
could arise or be already present. A recent clinical trial 
in India reported artesunate + SP efficacy dipping to near 
the WHO-established efficacy threshold of 90%, likely 
attributable to the presence of multiple pfdhfr and 
pfdhps mutations [53]. 

All samples harboured a single copy of the pfmdrl 
gene, suggesting that artesunate, lumefantrine, and mef- 
loquine are likely to have high efficacy in Pakistan, and 
that the latter two could be considered as alternate ACT 
partner drugs if the prevalence of SP resistance markers 
rises. Presence of point mutations at position N86Y, 
Y184F and N1042D of pfmdrl indicate that continued 
molecular monitoring of this gene may provide inform- 
ative insights on its association with artemisinin and 
partner drugs in this region. 

Conclusions 

The fixation of chloroquine resistance marker pfcrt 
K76T in Pakistan suggests that, in regions with known 
high levels of P, falciparum, improved species diagnosis 
and appropriate treatment with ACT are critical to en- 
sure that malaria patients receive effective care. The 
presence of multiple resistance alleles in pfdhfr and 
pfdhps demonstrates the utility of molecular surveillance 
to monitor the efficacy of artesunate + SP in Pakistan. 
Based on these findings, molecular and therapeutic 



efficacy studies to assess the continued efficacy of 
artesunate + SP are warranted. 

Competing interests 

The authors declare that they have r^o competir^g interests. 
Authors' contributions 

AAK designed study, carried out the laboratory experiments, and performed 
data analysis. AAK and MV drafted the manuscript. CGJ, MA and EMA 
contributed to laboratory experiments and data analysis. MFN and FN 
participated in sample data collection and microscopy. FH and TM 
contributed to sample collection and contributed advice for microsatellite 
data analysis. MV, SAM and CVP provided guidance and coordination for 
study design, genotyping, and data analyses, and edited and revised the 
manuscript. All authors read and approved the final manuscript. 

Acknowledgements 

The authors thank the study participants for their involvement in the study. 
The authors also acknowledge Matthew Adams, Malathi Vadia and 
Sudhaunshu Joshi at the University of Maryland for their technical assistance 
and Muhammad Shahid, Ismail Jalil, Ifthikhar Ahmed, Adeel Riaz and Noor 
Rehman in Pakistan for help in sample collection. The authors would like to 
thank Venkatachalam Udhayakumar at the Centers for Disease Control and 
Prevention for generously providing control DNA from Pakistan for 
microsatellite validation. This work was supported by a grant from the 
Higher Education Commission of Pakistan in support of AAK's PhD studies at 
Quaid-i-Azam University, Islamabad, Pakistan, and was supported in part by 
the Howard Hughes Medical Institute at the University of Maryland School of 
Medicine, Baltimore, MD, USA. 

Author details 

^Department of Biochemistry, Faculty of Biological Sciences, Quaid-i-Azam 
University, Islamabad, Pakistan. ^Howard Hughes Medical Institute/Center for 
Vaccine Development, University of Maryland School of Medicine, Baltimore, 
MD, USA. ^Worldwide Antimalarial Resistance Network Molecular Module, 
University of Maryland School of Medicine, Baltimore, MD, USA. ^King 
Edward Medical University, Lahore, Pakistan. ^Department of Biochemistry 
and Molecular Biology, University of Gujrat, Gujrat, Pakistan. ^Islamic 
International Medical College, Rawalpindi, Pakistan. ^Centre for Health 
Research, Divine Word University, Madang, Papua New Guinea. ^Department 
of Molecular and Cellular Parasitology, Juntendo University, Tokyo, Japan. 

Received: 30 May 2013 Accepted: 26 August 2013 
Published: 29 August 2013 

References 

1. Ghanchi NK, Ursing J, Beg MA, Veiga Ml, Jafri S, Martensson A: Prevalence 
of resistance associated polymorphisms in Plasmodium falciparum field 
isolates from southern Pakistan. Malar J 201 1, 10:18. 

2. Shah I, Rowland M, Mehmood P, Mujahid C, Razique F, Hewitt S, Durrani N: 
Chloroquine resistance in Pakistan and the upsurge of falciparum 
malaria in Pakistani and Afghan refugee populations, Ann Trop Med 
Paras/To/ 1997, 91:591-602. 

3. Bouma MJ, Dye C, van der Kaay HJ: Falciparum malaria and climate 
change in the northwest frontier province of Pakistan. Ann J Trop Med 
Hyg 1996, 55:131-137. 

4. Rab MA, Freeman TW, Durrani N, de Poerck D, Rowland MW: Resistance of 
Plasmodium falciparum malaria to chloroquine is widespread in eastern 
Afghanistan. Ann Trop Med Parasitol 2001, 95:41-46. 

5. Durrani AB, Durrani lU, Abbas N, Jabeen M: Epidemiology of cerebral 
malaria and its mortality. J Pak Med Assoc 1997, 47:213-215. 

6. WHO: World malaria report 201 1 Geneva: World Health Organization; 201 1 . 

7. Rowland M, Durrani N, Hewitt S, Sondorp E: Resistance of falciparum 
malaria to chloroquine and sulfadoxine-pyrimethamine in Afghan 
refugee settlements in western Pakistan: surveys by the general health 
services using a simplified in vivo test. Trop Med Int Health 1997, 
2:1049-1056. 

8. Howard N, Durrani N, Sanda S, Beshir K, Hallett R, Rowland M: Clinical trial 
of extended-dose chloroquine for treatment of resistant falciparum 
malaria among Afghan refugees in Pakistan. Malar J 20]], 10:171. 

9. WHO: World malaria report 2008. Geneva: World Health Organization; 2008. 



Khattak et at. Malaria Journal 201 3, 12:300 
http://www.nnalariajournal.conn/content/1 2/1 /300 



10. WHO: Division of communicoble disease control newsletter. Geneva: World 
Health Organization; 2005. 

1 1 . Fidock DA, Nomura T, Talley AK, Cooper RA, Dzekunov SM, Ferdig MT, Ursos 
LM, Sidhu AB, Naude B, Deitsch KW, Su XZ, Wootton JC, Roepe PD, Wellems 
TE: Mutations in the P. falciparum digestive vacuole transmembrane 
protein Pfcrt and evidence for their role in chloroquine resistance. 

Mol Cell 2000, 6:861-871. 

1 2. Wellems TE, Plowe CV: Chloroquine-resistant malaria. J Infect Dis 2001 , 
184:770-776. 

13. Djimde A, Doumbo OK, Cortese JF, Kayentao K, Doumbo S, Diourte Y, 
Coulibaly D, Dicko A, Su XZ, Nomura T, Fidock DA, Wellems TE, Plowe CV: A 
molecular marker for chloroquine-resistant falciparum malaria. N Engl J 
A/led 2001, 344:257-263. 

14. Wootton JC, Feng X, Ferdig MT, Cooper RA, Mu J, Baruch Dl, Magill AJ, Su 
XZ: Genetic diversity and chloroquine selective sweeps in Plasmodium 
falciparum. Nature 2002, 418:320-323. 

15. Rawasia WF, Sridaran S, Patel JC, Abdallah J, Ghanchi NK, Barnwell JW, 
Escalante AA, Udhayakumar V, Beg MA: Genetic backgrounds of the 
Plasmodium falciparum chloroquine resistant transporter {pfcrt) alleles in 
Pakistan. Infect Genet Evol 2012, 12:278-281. 

16. Mixson-Hayden T, Jain V, McCollum AM, Poe A, Nagpal AC, Dash AP, Stiles 
JK, Udhayakumar V, Singh N: Evidence of selective sweeps in genes 
conferring resistance to chloroquine and pyrimethamine in Plasmodium 
falciparum isolates in India. Antinnicrob Agents Chennother 2010, 
54:997-1006. 

17. DaRe JT, Mehlotra RK, Michon P, Mueller I, Reeder J, Sharma YD, Stoneking 
M, Zimmerman PA: Microsatellite polymorphism within pfcrt provides 
evidence of continuing evolution of chloroquine-resistant alleles in 
Papua New Guinea. Malar J 2007, 6:34. 

18. Cowman AF, Morry MJ, Biggs BA, Cross GA, Foote SJ: Amino acid changes 
linked to pyrimethamine resistance in the dihydrofolate reductase- 
thymidylate synthase gene of Plasmodium falciparum. Proc Natl Acad Sci 
USA 1988, 85:9109-9113. 

19. Peterson DS, Walliker D, Wellems TE: Evidence that a point mutation in 
dihydrofolate reductase-thymidylate synthase confers resistance to 
pyrimethamine in falciparum malaria. Proc Natl Acad Sci USA 1988, 
85:9114-9118. 

20. Peterson DS, Milhous WK, Wellems TE: Molecular basis of differential 
resistance to cycloguanil and pyrimethamine in Plasmodium falciparum 
malaria. Proc Natl Acad Sci USA] 990, 87:301 8-3022. 

21. Plowe CV, Djimde A, Bouare M, Doumbo 0, Wellems TE: Pyrimethamine 
and proguanil resistance-conferring mutations in Plasmodium falciparum 
dihydrofolate reductase: polymerase chain reaction methods for 
surveillance in Africa. Ann J Trop Med Hyg 1995, 52:565-568. 

22. Plowe CV, Kublin JG, Doumbo OK: P. falciparum dihydrofolate reductase 
and dihydropteroate synthase mutations: epidemiology and role in 
clinical resistance to antifolates. Drug Resist Updat 1998, 1:389-396. 

23. Kublin JG, Dzinjalamala FK, Kamwendo DD, Malkin EM, Cortese JF, Martino 
LM, Mukadam RA, Rogerson SJ, Lescano AG, Molyneux ME, Winstanley PA, 
Chimpeni P, Taylor TE, Plowe CV: Molecular markers for failure of 
sulfadoxine-pyrimethamine and chlorproguanil-dapsone treatment of 
Plasmodium falciparum malaria. J Infect Dis 2002, 185:380-388. 

24. Brooks DR, Wang P, Read M, Watkins WM, Sims PF, Hyde JE: Sequence 
variation of the hydroxymethyldihydropterin pyrophosphokinase: 
dihydropteroate synthase gene in lines of the human malaria parasite, 
Plasmodium falciparum, with differing resistance to sulfadoxine. 

EurJ Biochem 1994, 224:397-405. 

25. Foote SJ, Galatis D, Cowman AF: Amino acids in the dihydrofolate 
reductase-thymidylate synthase gene of Plasmodium falciparum involved 
in cycloguanil resistance differ from those involved in pyrimethamine 
resistance. Proc Natl Acad Sci USA 1 990, 87:3014-301 7. 

26. Triglia T, Cowman AF: Primary structure and expression of the 
dihydropteroate synthetase gene of Plasmodium falciparum. Proc Natl 
Acad Sci USA 1 994, 91 :71 49-71 53. 

27. Triglia T, Menting JG, Wilson C, Cowman AF: Mutations in dihydropteroate 
synthase are responsible for sulfone and sulfonamide resistance in 
Plasmodium falciparum. Proc Natl Acad Sci USA 1 997, 94:1 3944-1 3949. 

28. Wang P, Read M, Sims PF, Hyde JE: Sulfadoxine resistance in the human 
malaria parasite Plasmodium falciparum is determined by mutations in 
dihydropteroate synthetase and an additional factor associated with 
folate utilization. Mol Microbiol 1997, 23:979-986. 



Page 8 of 9 



29. Happi CT, Gbotosho GO, Folarin OA, Akinboye DO, Yusuf BO, Ebong 00, 
Sowunmi A, Kyle DE, Milhous W, Wirth DF, Oduola AM: Polymorphisms in 
Plasmodium falciparum dhfr and dhps genes and age related in vivo 
sulfadoxine-pyrimethamine resistance in malaria-infected patients from 
Nigeria. Acta Trop 2005, 95:183-193. 

30. Picot S, Olliaro P, de Monbrison F, Bienvenu AL, Price RN, Ringwald P: A 
systematic review and meta-analysis of evidence for correlation 
between molecular markers of parasite resistance and treatment 
outcome in falciparum malaria. Malar J 2009, 8:89. 

31. Staedke SG, Kamya MR, Dorsey G, Gasasira A, Ndeezi G, Charlebois ED, 
Rosenthal PJ: Amodiaquine, sulfadoxine/pyrimethamine, and 
combination therapy for treatment of uncomplicated falciparum malaria 
in Kampala, Uganda: a randomised trial. Lancet 2001, 358:368-374. 

32. Lim P, Wongsrichanalai C, Chim P, Khim N, Kim S, Chy S, Sem R, Nhem S, Yi P, 
Duong S, Bouth DM, Genton B, Beck HP, Gobert JG, Rogers WO, Coppee JY, 
FandeurT, Mercereau-Puijalon 0, Ringwald P, Le BJ, Ariey F: Decreased in vitro 
susceptibility of Plasmodium falciparum isolates to artesunate, mefloquine, 
chloroquine, and quinine in Cambodia from 2001 to 2007. Antinnicrob 
Agents Chemother 2010, 54:2135-2142. 

33. Price RN, Uhlemann AC, Brockman A, McGready R, Ashley E, Phaipun L, 
Patel R, Laing K, Looareesuwan S, White NJ, Nosten F, Krishna S: Mefloquine 
resistance in Plasmodium falciparum and increased pfmdri gene copy 
number. Lancet 2004, 364:438-447. 

34. Price RN, Uhlemann AC, van VM, Brockman A, Hutagalung R, Nair S, Nash D, 
Singhasivanon P, Anderson TJ, Krishna S, White NJ, Nosten F: Molecular and 
pharmacological determinants of the therapeutic response to 
artemether-lumefantrine in multidrug-resistant Plasmodium falciparum 
malaria. Clin Infect Dis 2006, 42:1570-1577. 

35. Alker AP, Lim P, Sem R, Shah NK, Yi P, Bouth DM, Tsuyuoka R, Maguire JD, 
FandeurT, Ariey F, Wongsrichanalai C, Meshnick SR: Pfmdri and in vivo 
resistance to artesunate-mefloquine in falciparum malaria on the 
Cambodian-Thai border. Am J Trop Med Hyg 2007, 76:641-647. 

36. Petersen I, Eastman R, Lanzer M: Drug-resistant malaria: molecular 
mechanisms and implications for public health. FEBS Lett 2011, 
585:1551-1562. 

37. Khatoon L, Baliraine FN, Bonizzoni M, Malik SA, Yan G: Prevalence of 
antimalarial drug resistance mutations in Plasmodium vivax and P. 
falciparum from a malaria-endemic area of Pakistan. Am J Trop Med 
Hyg 2009, 81:525-528. 

38. Zakeri S, Gil JP, Bereckzy S, Djadid ND, Bjorkman A: High prevalence of 
double Plasmodium falciparum dhfr mutations at codons 108 and 59 in 
the Sistan-Baluchistan province, Iran. J Infect Dis 2003, 187:1828-1829. 

39. Nawaz S: EAJA—A Most Dangerous Place. Washington, D.C; 2009. 

40. Asif SA: Departmental audit of malaria control programme 2001-2005 
north west frontier province {mJV?).J Ayub Med Coll Abbottabad 2008, 
20:98-102. 

41. Snounou G, Viriyakosol S, Jarra W, Thaithong S, Brown KN: Identification of 
the four human malaria parasite species in field samples by the 
polymerase chain reaction and detection of a high prevalence of mixed 
infections. Mol Biochem Parasitol 1993, 58:283-292. 

42. Zhou Z, Poe AC, Limor J, Grady KK, Goldman I, McCollum AM, Escalante AA, 
Barnwell JW, Udhayakumar V: Pyrosequencing, a high-throughput method 
for detecting single nucleotide polymorphisms in the dihydrofolate 
reductase and dihydropteroate synthetase genes of Plasmodium 
falciparum. J Clin Microbiol 2006, 44:3900-3910. 

43. CVD Malaria Group: Protocols, [http://medschool.umaryland.edu/malaria/ 
protocols.asp]. 

44. Nash D, Nair S, Mayxay M, Newton PN, Guthmann JP, Nosten F, Anderson 
TJ: Selection strength and hitchhiking around two anti-malarial 
resistance genes. Proc Biol Sci 2005, 272:1 153-1161. 

45. Laufer MK, Takala-Harrison S, Dzinjalamala FK, Stine OC, Taylor TE, Plowe CV: 
Return of chloroquine-susceptible falciparum malaria in Malawi was a 
reexpansion of diverse susceptible parasites. J Infect Dis 2010, 202:801-808. 

46. Lumb V, Madan R, Das MK, Rawat V, Dev V, Khan W, Khan H, Sharma 
YD: Differential genetic hitchhiking around mutant pfcrt alleles in the 
Indian Plasmodium falciparum population. J Antimicrob Chemother 
2012, 67:600-608. 

47. Ahmed A, Bararia D, Vinayak S, Yameen M, Biswas S, Dev V, Kumar A, Ansari 
MA, Sharma YD: Plasmodium falciparum isolates in India exhibit a 
progressive increase in mutations associated with sulfadoxine- 
pyrimethamine resistance. Antimicrob Agents Chemother 2004, 48:879-889. 



Khattak et at. Malaria Journal 201 3, 12:300 
http://www.nnalariajournal.conn/content/1 2/1 /300 



Page 9 of 9 



48. Biswas S, Escalante A, Chaiyaroj S, Angkasekwinai P, Lai AA: Prevalence of 
point mutations in the dihydrofolate reductase and dihydropteroate 
synthetase genes of Plasmodium falciparum isolates from India and 
Thailand: a molecular epidemiologic study. Trop Med Int Health 2000, 
5:737-743. 

49. Zakeri S, Farahani MS, Afsharpad M, Salehi M, Raeisi A, Djadid ND: High 
prevalence of the 437G mutation associated with sulfadoxine resistance 
among Plasmodium falciparum clinical isolates from Iran, three years 
after the introduction of sulfadoxine-pyrimethamine. Int J Infect Dis 2010, 
14(Suppl 3):e123-el28. 

50. Hapuarachchi HC, Dayanath MY, Bandara KB, Abeysundara S, 
Abeyewickreme W, de Silva NR, Hunt SY, Sibley CH: Point mutations in the 
dihydrofolate reductase and dihydropteroate synthase genes of 
Plasmodium falciparum and resistance to sulfadoxine-pyrimethamine in 
Sri Lanka. Am J Trop Med Hyg 2006, 74:198-204. 

51. Parikh R, Amole I, Tarpley M, Gbadero D, Davidson M, Vermund SH: Cost 
comparison of microscopy vs. empiric treatment for malaria in 
Southwestern Nigeria: a prospective study. Malar J 20]0, 9:371. 

52. Kakar Q, Khan MA, Bile KM: Malaria control in Pakistan: new tools at hand 
but challenging epidemiological realities. East Mediterr Health J 2010, 
16(Suppl):S54-S60. 

53. Saha P, Guha SK, Das S, Mullick S, Ganguly S, Biswas A, Bera DK, 
Chattopadhyay G, Das M, Kundu PK, Ray K, Maji AK: Comparative efficacies 
of artemisinin combination therapies in Plasmodium falciparum malaria 
and polymorphism of pfATPased, pfcrt, pfdhfr, and pfdhps genes in tea 
gardens of Jalpaiguri District, India. Antimicrob Agents Chemother 2012, 
56:2511-2517. 



doi:1 0.1 1 86/1 475-2875-1 2-300 

Cite this article as: Khattak et al.: A comprehensive survey of 
polymorphisms conferring anti-malarial resistance in Plasmodium 
falciparum across Pakistan. Malaria Journal 2013 12:300. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at Ri^nHod Tpntrai 

www.biomedcentral.com/submit momea L,enirai 



